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Fluid Terpperature Measurement 



The present invention relates to a device and method for measuring the temperature of a 
fluid. It has particular application for total temperature measurement in high temperature 
5 unsteady gas flows, for example in gas tuVbine engines or similar high temperature 
environments. 

Accurate temperature measurement within the aggressive environment of a gas turbine 
engine presents significant problems. Combustor exit temperatures typically reach 

10 around 2000^K and will rapidly^destroy any conventional form of instrumentation exposed 
to that environment. Thermocouples are used in many applications to provide low cost 
measurement of high temperatures, but their time response is poor and restricts their 
usefulness for measurement within unsteady flows. Neither could they survive in a 
typical gas turbine combustor exit for long enough to provide reliable temperature 

IS measurements. 



GB-A-2314164 describes an alternative form of device for such measurement comprising 
a probe with two thin film heat transfer gauges at its tip. Each gauge comprises a 
platinum element whose electrical resistance varies with temperature and is sensed 
20 through an associated circuit Prior to exposure to the environment to be measured one 
of the gauges is preheated, by use of an external heater, a separate heating element 
incorporated in the probe, or by passing a specified current through the respective 
platinum element. In use the probe is rapidly inserted into and retracted from the 
environment and the fluid temperature is determined from the relationship: 

25 

Tt == Twi + qi(Tw2 - Twi) / (qi - qa) 

where T| is the total temperature of the.fluld being measured, Twi and T^a are the 
respective suri'ace temperatures of the two gauge elements arising from exposure to the 

30 fluid (and consequent on their initial temperature difference) and qi and o^. are the 
respective heat transfer rates experienced by the two gauge elements (which can be 
derived from their temperature histories). The requirement to preheat one of the 
elements adds undesirable complexity to the process, however. It also places a practical 
limitation on the temperatures which can be measured, as the initial differential is rapidly 

35 eliminated as the two elements are heated by the fluid. 
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The present invention therefore seelcs to provide a means of temperature measurement 
for the kind of service Indicated above which avoids the draWbacl<s of preheating as in 
GB-A-2314164. 

5 In a first aspect the invention resides in a device for use in measuring the temperature of 
a fluid comprising a structure bearing two temperature sensitive elements adapted to be* 
temporarily exposed to the fluid, wherein the structure provides respective regions for the 
diffusion of heat from the fluid through the respective said elements, the themnal products 
of said regions being selected such that, in use. said elements experience different heat 
10 transfer rates when exposed to the same fluid temperature. 

In a second aspect the invention resides in a method of measuring the temperature of a 
fluid which comprises temporarily exposing to the fluid a device in accordance with the ' 
first aspect of the invention; monitoring the respective temperatures of the temperature 
15 sensitive elements of such device over a period; deriving from respective changes of 
temperature of said elements the respective heat transfer rates experienced thereby; and 
deriving the temperature of the fluid from a relationship of the temperatures of said 
elements and the derived heat transfer rates. 

20 In a third aspect the invention resides in an apparatus for measuring the temperature of a 
fluid comprising a device according to the first aspect of the invention; means for 
monitoring the respective temperatures of the temperature sensitive elements of such 
device over a period; and computational means for deriving from respective changes of 
temperature of said elements the respective heat transfer rates experienced thereby and 

25 for deriving the temperature of the fluid from a relationship of the temperatures of said 
elements and the derived heat transfer rates. 

in use OT ine presem invention nign Tiuia temperatures can be derived from 
measurements of the surface temperatures of the structure of the exposed device at the 

30 locations of the two temperature sensitive elements, together with the derived associated 
heat transfer rates, without the need for preheating one of the elements, since the nature 
of the structure itself will give rise to a differential between tiie two locations. In addition 
to simplifying the measurement process this removes the limit on measurable fluid 
temperatures inherent due to preheating in GB-A-2314164. the only practical limitation 

35 being the survivability of the materials used in tiie device. 
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These and other features of the invention will now be more particularly described, by way 
of exannple only, with reference to the accompanying drawings, in which: 

Figure 1 is a side view of one preferred embodiment of a temperature measurement 
5 probe according to the invention; 

Figure 2 is a graph of surface temperatures on a probe as illustrated in Figure 1 as 
measured under test conditions, together with the fluid temperature as measured by a 
separate thermocouple; 

10 

Figure 3 is a graph of heat transfer rates as derived from the probe surface temperature 
measurements illustrated In Figure 2; 

Figure 4 is a graph of flow total temperature as derived from the probe surface 
15 temperature measurements illustrated in Figure 2, and as measured by the 
thermocouple; and 

Figures 5 and 6 are side views of two further embodiments of a temperature 
measurement probe according to the invention. 

20 

Figure 1 illustrates a probe P1 for temperature measurement in accordance with the 
invention. It comprises a hollow steel shaft 1 canying a circular cross-section rod 2 of 
electrically insulative material, in this example a machinable glass known under the 
registered trade mark Macor (of Corning Incorporated). The distal end of the rod 2 Is 

25 coaxially machined. to leave a thin cylindrical wall 3 of Macor surrounding an air pocket 4, 
and is closed by a Macor end cap 5. In one example which has been tested as 
described below |he diameter of the rod 2 js 4mm and the wall thickness at 3 is 0.25mm. 
Two substantially identical platinum thin film heat transfer elements (resistance 
thermometers) 6 and 7, typically about 1^m thick, are painted or sputtered in the 

30 circumferential direction onto the Macor surface near the end of the rod so that one of 
them (6) lies on the thin wail 3 while the other (7) lies on the full thickness of the rod 2. 
Electrical leads are connected to the ends of the respective elements 6 and 7 for the 
supply of a constant current to each. These leads are not shown in Figure 1 but may 
comprise tracks of conductive ceramic paste painted onto the rod 2 and extending to a 

35 connector within a Macor collar 8 at the junction with the steel shaft 1 , from where a 
conventional cable is taken through the shaft to an external power/monitoring circuit. 
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In use of the illustrated probe to measure the total temperature of a high temperature gas 
flow, such as within part of a gas turbine engine, it is temporarily inserted into the stream 
with the elements 6 and 7 facing the direction of flow (i.e. with the longitudinal axis of the 
5 probe transverse to the flow), and withdrawn again before the rod 2 reaches its softening 
temperature. The platinum elements 6 and 7 are powered by a respective supply of 
constant current and their temperatures over the period of exposure are measured by 
monitoring the resultant voltages across them. 

10 The elements 6 and 7 are located close to the stagnation point of the cylindrical probe 
surface. Therefore the convective heat transfer between the flow and the probe is 
proportional to the temperature difference between the flow total temperature and the 
probe surface temperature, such that: 

15 qi=hi(Tt-Twi) and q2 =F h2 (Tt - Tw2) 

where qi and q2 are the surface heat transfer rates (W/m^) at the respective locations of 
the elements 6 and 7, hi and h2 are the convective heat transfer coefficients (W/m^K) at 
those locations, Ti is the flow total or stagnation temperature and T^i and Tw2 are the 
20 probe wall temperatures at the respective locations. The convective heat transfer 

coefficient h is a property, of the stagnation point boundary layer and is thus a function of 
the probe geometry and the upstream flow condition. Since it can be assumed that at 
the locations of the elements 6 and 7 two geometrically identical cylinders are exposed to 
the same flow then hi=h2 and the foregoing equations can be resolved to: 

25 

Ti = Twi + qi (Tw2 Twi)/(qi-q2) 

In use, as indicated above, the surface temperatures Twi and Tw2 are measured by 
monitoring the voltages across the respective thin film elements 6 and 7. which typically 

30 have a frequency response up to lOOKHz. Tt can be calculated from the above 
relationship for successive instants during the period of exposure, in an associated 
computer, at intervals selected with regard to the frequency response of the thin films 
and the nature of the flow concerned - typically at 10ms intervals for gas turbine flow 
investigations. For this purpose the convective heat fluxes qi and qz at each instant are 

35 derived from the thin film temperature histories during the period of exposure up to that 
instant using a transient one-dimensional heat diffusion process which models the 
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transient heat diffusion within the probe structure underlying the respective elements, as 
will be more fully explained hereafter. 

Clearly, measurement of Tt by this method requires the establishment of different q 
5 values at the two Tw locations, which in the illustrated probe is achieved by the presence 
of different materials in the regions of the probe structure into which the heat diffuses 
from those locations - in the case of element 7 solid Macor and in the case of element 6 
a thin wall of Macor followed by air. In the probe design the important parameter is the 

thermal product ^pck of the materials used, where p is density, c is specific heat 

10 capacity and k is conductivity. The ./pcfc .of IWacor is 1 61 5 J/m^Ks^'^ whereas for air it 
is 5.5J/m^Ks'*'^. Upon initial exposure of the probe to the high temperature fluid the 
temperature response of both elements 6 and 7 will be the same as both are present on 
a Macor substrate. However, after a short period given by t = x^/a, where x is the 
thickness of the wall 3 and a^klpc - and typically 10ms for a Macor wall thiclcness of 

15 0.25mm - the element 6 will experience a themnal response related to the air in pocket 4. 
This results in a higher surface temperature and a lower heat transfer rate than for the 
element 7. ' As the surface temperatures increase vAXh time the heat transfer rates 
reduce, since the heat transfer rate is dependent on the temperature difference between 
the flow and the probe surface. 

20 

Turning now to the detemiination of the surface heat fluxes, it is known that for a film on 
a flat plate the following semi-infinite flat plate heat conduction equation can be used: 

d^T ^ 1 dT 
dx^ a dt 

25 

where T is the temperature, x is the dlistance measured from the surtace, t is the time 
and a=k/pci with the following boundary conditions applied: 

T(oo,t) = T|, initial temperature of substrate 
30 T(0,t) = Ts(t), measured surface temperature history. 

This heat conduction equation can be solved for surface heat transfer rate using either 
electric analogue circuits or numerically. Methods commonly used include Numerical 
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Integration (Sdiultz and Jones, 1973), Discrete Fourier Transform (Moss and Ainsworth 
1993), and the Impulse Response Method of Oldfield (2000). 

However, in the case of the illustrated probe the thin films 6 arid 7 are mounted on a 
5 cylindrical surface. For heat transfer gauges on geometries where the surface curvature 
is comparable to the heat penetration depth during the operation time, allowance has to 
be made for the curvature effect. This can either be based on correction factors for a flat 
plate solution or radius variation has to be included in the one-dimensional heat 
conduction equation. Correction methods only perform well for short time periods of run 
10 time and low heat transfer rates. Furthermore, in the illustrated probe significant 
variations are found in the substrate thermal properties. To take account of these 
variations they have to be incorporated in the one-dimensional heat conduction equation. 
The methods for solving the heat conduction equation given above do not allow for the 
radius and thermal property variation to be taken into account. They are developed 
15 specifically for flat plate solutions with fixed substrate properties. 

We have therefore derived the following heat conduction equation for both radius of 
curvature and thermal property variation: 



20 



dr 



dr 



r dr at 



where r Is the radius within the substrate, p the density, k conductivity and c specific heat 
capacity, with boundary conditions: 



25 — = 0 , i.e. symmetry applies at the centre of the substrate 

L9^J(o.o 

T(R,t) = Ts(t), measured surface temperature history 



where R Is the radius at the surface of the probe. 

30 

The symmetry boundary condition does not restrict tiie heat transfer solution to be 
obtained once the heat penetration has reached the centre of the substrate. In practice 
the solution is limited by the time taken for two and three-dimensional effects to become 
significant. These are associated witti the physical constiruction of the probe body: 



'35 
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Various methods can be used to solve this heat conduction equation for dTldt (which 
gives q), the preferred method being based on a finite difference scheme. To achieve 
this the probe is meshed from the centre to the surface with nodes in the radial direction. 
The spacial temperature within the probe body is derived in time using the measured 
surface temperature. The themial properties are evaluated at each time step for the 
previously evaluated temporal temperature. Since the solution requires the inversion of a 
matrix, it is desirable to minimise the number of nodes within the probe to reduce 
analysis time. Hence, grid refinement is employed at tiie surface of flie probe where ttie 
temperature gradients will be the largest The grid refinement also takes into account Oie 
diameter of the probe. In addition ttiis method allows the use of multi-layered substrates. 
For this case the appropriate thermal properties are calculated for each layer as the 
solution progresses though the probe substrates. The grid refinement in this case is 
optimised for both the substi^te materials. The substrates are not restricted to solids, 
they may include gases or gas mixtures (such as the air pocket 4) so long as their 
thermal properties are known, or even a vacuum. 

Figures 2 to 4 show test results for a probe P1 substantially as described with reference 
to Figure 1 which has been exposed substantially instantaneously to a fluctuating gas jet 
between about 500 and 600^K. Figure 2 shoWs the probe surface temperature rise as 
measured by tiie elements 6 and 7 and a trace from a fast response tiiermocouple 
placed adjacent to tiie probe in the same jet. Figure 3 shows the calculated heat 
transfer rate, from the surface temperature history, plotted against time for each element 
6.7. Figure 4 shows the evaluated flow total temperature from the tiiin film 
measurements of the Figure 1 probe and the total temperature as measured by the 
thermocouple. Note that the fluctuations on the probe trace are real as the jet used in 
the test was fluctuating. The thenmocduple does not, however, have suffident frequency 
response to captijre the fluctuations. . 

The illustrated probe measures the flow total temperature in-espective of the flow Mach 
number because stagnation entiialpy is conserved as Uie flow decelerates to tiie 
stagnation point and tiiere is virtually no viscous dissipation in the stagnation point 
boundary layer. Indeed, even in flows where an immersed body would eventually reach 
a recovery temperature slightiy lower tiian tiie flow total temperature, the probe 
accurately measures the total temperature. This is because the thin films 6 and 7 
measure the transient heat flux close to the stagnation point, which is driven by the flow 
total temperature. 
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In a variant of the probe P1 the rod 2 Is made from a refractory ceramic, such as 
giumina, having a substantially higher melting point than Macor, and hence capable of 
use in higher temperature flows and/or over a longer exposure time. Another candidate 
5 material is quartz. 

Figures 5 and 6 illustrate alternative probe constmctions P2 and P3, with like reference 
numerals denoting like parts with Rgure 1 . in Figure 5 the rod 2 is replaced with a tube 
9 of IVIacor, ceramic or quartz sun-ounding a metal rod 10 over most of its length but still 

10 leaving an air pocket 4 at its distal end: The rod 10 has a significantly higher thenmal 
product than the tube 9. In this case the thin film element 7 lies on a portion of the tube 
9 overlying the rod 10, which increases the heat transfer rate into the structure and leads 
to a greater temperature differential between the two thin films than in the Figure 1 
embodiment. The rod 10 also enhances the physical strength of the probe. In Figure 6 

15 the rod 2 is replaced with a solid rod composed of two axial sections 1 1 and 12 made of 
respective materials with different thermal products (e.g. Macor and ceramic or quartz 
and ceramic), and one of the thin film elements 6,7. is applied to each section. 

Altiiough in each embodiment described above the thin film elements 6,7 are shown 
20 deposited on the respective substrate in its circumferential direction, they could 
alternatively be deposited in the longitudinal direction. 



